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Introduction générale 


Introduction générale. 

Les structures dites « polymères cellulaires » occupent une place importante dans de 
nombreux domaines industriels. Elles sont présentes pratiquement dans tous les domaines 
d’application ; ces structures résultent de la volonté d’introduire une certaines proportion de 
cavités dans un squelette de polymère dans le but : 

-d’augmenter la rigidité et la légèreté (utilisation divers). 

-d’améliorer le caractère isolant, thermique ou phonique (mousses pour bâtiment), 
-d’augmenter la capacité d’absorber de l’énergie en compression (mousse d’emballage). 

La théorie des déformations élastiques appliquée aux structures isotropes n’est plus valable 
dans le cas des solides cellulaires, ceci est dû à la complexité de leur comportement 
mécanique grâce à leur morphologie. Deux voies principales qui permettent de traiter 
partiellement leur comportement sont, actuellement, largement utilisées : 

La première se situe dans une approche microscopique, repose sur l’analyse dimensionnelle. 

Plusieurs travaux ont été effectués concernant l’analyse dimensionnelle des matériaux 
cellulaires : 

Ko [1] , Gent et Thomas [2 ' 3 4] , Gent et Rusch [5 ’ 6] , Patel et finnie [15] , Dementjev et Tarakanov 
[7] , Zhu et al [8 9] , Hilyard [10] , Gibson et Ashby [30] , Zhang et Ashby [11] , Ould Eyih [12] , Pradel 
[13] . La majorité de ces études présente les mousses comme un assemblage de poutres 
représentant les arêtes des cellules, avec des faces dans le cas des mousses à cellules fermées, 
en utilisant des structures géométriques différentes telles que les structures cubiques, 
tétraédrique, dodécaèdrique, tétrakaidécaèdrique. 

La seconde voie consiste à l’analyse macroscopique sous une approche numérique par la 
méthode des éléments finis. 

Dans ce contexte on cite les études suivantes : Renz [1617] et Czysz [18] Où ces derniers auteurs 
ont décrit la phase élastique par la loi de Hook. 

Il est à noter que ces auteures se sont intéressées à la partie linéaire et ont négligé la partie non 
linéaire. 
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Actuellement, les modèles non linéaires pour les solides cellulaires sont, en général, fondés 
sur une approche hyperélastique pour le cas compressible avec une énergie de déformation 
particulière. 

La plupart de ces études sont fondées sur une fonction énergie similaire à celle utilisée dans 
le modèle d’Ogden [19 ], à la différence qu’elle est adaptée au cas d’un matériau compressible. 
Cette fonction a été proposée par Hill |2X| , storakers [29] et ogden [19] et introduite par plusieurs 
auteurs dans l’outil- FE ABAQUS afin de décrire les propriétés mécaniques des solides 
cellulaires dans le cas de la compression uniaxiale. 

Le but de ce présent travail consiste en l’étude comportementale mécanique des polymères 
cellulaires par la méthode des éléments finis. La mousse étudiée est le polyuréthane à cellules 
ouvertes. 

Le premier chapitre de ce mémoire expose les polymères cellulaire, dites aussi alvéolaires 
leur nature, leurs multiples propriétés, leurs diversités d’applications ainsi que les bases de 
chimie de la mousse de polyuréthane et les différents procédés de fabrication. 

Le chapitre deux résume une analyse dimensionnelle des solides cellulaires de structures 
périodiques de différentes géométries. 

Le troisième chapitre est dédié à la simulation du comportement en compression uniaxiale qui 
nous permettra d’étudier et visualiser la réponse mécanique des mousses de polyuréthane. Les 
résultats obtenus représentent dans leur majorité une comparaison entre la mousse et le 
matériau massif qui la compose afin de mettre en évidence la capacité des mousses d’absorber 
de l’énergie en compression. L’étude de l’influence de la densité sur le début de la 
densification est à envisager. Le cas d’une structure trouée a permis de bien visualiser 
l’incidence des ouvertures sur la rigidité du matériau. 
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Chapitre 1 


Les polymères cellulaires. 


1.1. Introduction. 

Les mousses polymères dites aussi alvéolaires sont constitués d’un squelette de polymère 
entourant des alvéoles [20] . Elles sont présentes pratiquement dans tous les domaines 
d’application. Elles peuvent être produites à partir d’une grande variété de polymères 
synthétiques, dont le polychlorure de vinyle (PVC) les polyoléfines (PE, PP), le polystyrène 
(PS), le polyuréthane (PUR), le polyméthacrylimide (acrylique, PMI), les phénoliques (PF), 
le polyéthérimide (PEI) ou les copolymères styrène-acrylonitrile (SAN). Elles sont 
disponibles dans une gamme de masses volumiques allant de moins de 30 Kg/m à plus de 
300 Kg/m 3 [20 ’ 44] . 

1.2. Types des structures cellulaires. 

Nous pouvons décrire les solides cellulaires comme un assemblage de cellules vides, chacune 
étant entourée par des arêtes tels que les solides cellulaires à porosité ouverte (figure 1.1.a) ou 
des parois tels que les solides cellulaires à porosité fermée (figure 1.1.b) [41] . 


Coupe transversale de l'arête 

Arête 

Sommet 



Coupe transversale 


Coupe transversale 


Fig. 1.1 : Image par MEB de (a) polyuréthane à cellules ouvertes ; (b) a cellules fermées [37] . 
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1.3. Fonctionnalités des polymères cellulaires [20] . 

Dans les polymères cellulaires, le taux de gaz plus ou moins piégé dans les alvéoles à parois 
très minces entraine des propriétés physiques, mécaniques et chimiques différentes de celles 
des mêmes polymères compacts ; ces propriétés des solides cellulaires ouvrent un grand 
nombre d’applications émergentes tels que : 

-amélioration du caractère isolant, thermique ou phonique : mousse pour le bâtiment ; 

-amélioration du pouvoir amortissant : mousses d’emballage ; 

-allègement avec conservation de certaines propriétés structurelles ; 

-augmentation de la souplesse : premières mousse de latex pour l’ameublement ; 

-diminution de la masse volumique permettant d’occuper des espaces avec des masses de 
matériau réduites ; 

-dégradation de la résistance au vieillissement, de la tenue chimique et plus grande sensibilité 
au feu à cause de la division de la matière ; 

-dégradation des valeurs absolues des propriétés mécaniques même si les valeurs ramenées à 
la masse sont intéressantes. 

Parmi les fonctionnalités les plus courantes, citons : 

-l’amortissement des chocs et vibration ; 

-l’isolation thermique et phonique ; 

-l’étanchéité ; 

-la flottabilité ; 

-le remplissage d’espace avec faibles masses ; 

-les âmes de composites ; 

-l’allègement des matériaux structuraux ; 

-la filtration, l’absorption. 
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1.4. Propriétés des polymères cellulaires [41] . 

1.4.1. Propriétés électriques. 

La résistivité électrique d’une mousse R* fait avec un matériau conducteur de résistivité R s est 
inversement proportionnelle à la densité. Si l’on considère une mousse à cellules cubiques 
fermées, la densité relative sera de : 


Rs__ 2 
R* ~ 3 PS 


Pour une mousse à cellule ouverte : 


Rs__ 1 £_ 
R* 3 Ps 


( 1 . 1 ) 


( 1 . 2 ) 


La résistivité d’un solide cellulaire ne dépend pas seulement de la densité mais aussi de la 
géométrie des pores. 

1.4.2. L’isolation thermique. 

En général il ya trois modes de transport de la chaleur : la convection, la conduction et la 
radiation. Dans les mousses à cellules fermées que l’on utilise pour l’isolation thermique, la 
convection ne joue généralement aucun rôle. La conduction thermique est principalement due 
aux faces et des arêtes de la structure, mais le gaz enfermé dans les pores peut y contribuer de 
façon non négligeable. On peut donc estimer que le coefficient de conduction thermique sera 

une fonction linéaire de la densité relative de la mousse, — . Pour estimer le transfert de 

PS 

chaleur par radiation, on peut partir de l’idée que l’énergie émise sur une face de la couche 
sera réabsorbée par l’autre. L’énergie émise est donnée par la loi de Stephan-Boltzmann, ce 
qui permet d’écrire le flux d’énergie : 

<7 = <i 0 UT + AT)*-T*) = 4 </ 0 r J Ar (1.3) 


Où : q 0 est une constante de Stephan- Boltzmann. Mais les rayons infrarouge doivent aussi 
passer à travers les parois du matériau cellulaire, une partie sera absorbée, ce qui signifie que 
l’équation (1.3) doit être multipliée par un facteur de transmission de forme : 

exp { -KBp/p^ ) ( 1 -4) 
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K : étant une constante décrivent l’absorption dans les parois cellulaires. 

Finalement, tenant compte de la conduction ainsi que du transport par radiation, on peut écrire 
le coefficient de transport thermique de la couche isolante : 

A ' = + 4< yo^ 3 *cxp|- A ^ | (1.5) 


Les deux premiers termes décrivent le transport par conduction, respectivement dans les 
parois et dans le gaz remplissant les cellules, et le troisième décrit le transport par radiation. 

{ A s \ la conductivité du gaz dans les pores. 

A g \ la conductivité du matériau massif. 

B\ l'épaisseur d'une couche isolante. 

T\ la température en kalvin. 

1.4.3. L’absorption acoustique. 

Les matériaux cellulaires sont inutilisables pour cette tache : lorsque les mousses absorbent 
les ondes sonores, le bruit peut être transporté par le mur lui-même, par exemple à cause 
d’une source de vibration agissant directement sur le mur ; ce genre de vibration peut bien 
être absorbé par les matériaux cellulaires surtout à cellules ouvertes. 

1.4.4. L’absorption de l’énergie en compression. 


Les solides cellulaires, notamment les alvéolaires, sont capable d’absorber de l’énergie en 
compression à cause du plateau dans les courbes contrainte-déformation en compression. En, 
effet, ce plateau qui est dû au fléchissement réversible des parois, permet d’absorber une large 
partie de l’énergie de déformation. 

1.4.4.1. Essai de compression /décompression. 

Pour les matériaux alvéolaires tels que les mousses de polyuréthane, les propriétés 
mécaniques sont souvent caractérisées par des tests de compression suivis d’une 
décompression [ ] . Les échantillons sont couramment comprimés entre les deux plateaux de 
compression de l’appareil jusqu'à une déformation souhaitée (figure 1.2) L’échantillon subit 
ensuite une décompression jusqu'au retour à la positon initiale. 


6 








Chapitre 1 


Les polymères cellulaires. 


Contrainte 


Plaques de 
compression 



Fig. 1.2 : Essai de compression /décompression [56] . 

1.5. Le polyuréthane cellulaire. 

Les polyuréthanes, polymères thermodurcissables occupent une place privilégiée dans 
l’industrie [46] . 

En 1937 le Dr Otto Bayer et ses collaborateurs ont réalisé la synthèse du premier 
polyuréthane. Ils ont fait réagir directement le butanediol sur l’hexaméthylène isocyanates 
(figure 1.3). 


HO 


OH 


OC N 


NI CO 



Fig. 1.3 : Polyaddition [54] . 

Depuis leur découverte, les polyuréthanes sont imposés dans tous les domaines de notre vie 
quotidienne grâce à leur nature, leurs multiples propriétés, notamment de confort et 
d’isolation, et leur diversité d’applications [45] . Les polyuréthanes sont des polymères 
contenant des groupes 


7 






Chapitre 1 


Les polymères cellulaires. 


uréthanes (appelés également carbamates) dans la chaîne moléculaire (figure 1.4). Outre le 
groupe uréthane, ces matériaux contiennent aussi plusieurs autres types de groupes, tels que 
l’amide, l’urée, l’éther, et l’ester [42] ’ [43] . 

-fO-R-O-C-N-R'-N-C-Vr 

II I I II 

OH HO 

Fig. 1.4 : Formule d’un polyuréthane linéaire [54] . 

1.5.1. Applications de polyuréthanes cellulaires. 

Grâce à leur basse densité, et ce que l’on voit bien dans la figure (1.5) les mousses de 
polyuréthanes occupent une place importante dans tous les domaines de construction tels que 
l’emballage, l’ameublement et le transport [42] . Ces mousses présentent des températures 
d’utilisation relativement élevées (150°C), une bonne capacité d’absorption acoustique, et 
peuvent être aisément découpées et usinées aux géométries requises [47] . 

Les domaines d’application des mousses polyuréthanes sont nombreux [46] et cette liste n’est 
pas exhaustive : 

-colles : (résistante à l’eau) : la mousse de polyuréthane est à titre un excellent matériau de 
collage et d’isolation thermique utilisés par exemple dans le bâtiment (collage de panneaux 
d’isolation sur des cloisons) ; 

-amortisseurs : le polyuréthane chargés de composantes anti- abrasion (roue et roulettes pour 
patins ou planche à roulettes, caddy, etc.), polyuréthane micro alvéolaires (silentblocs et par 
chocs) ; 

-literie, matelasserie ; 

-décoration : fausses poutres, sculptures,.qui sont ensuite peintes ; 

- nautisme : de par ses qualités d’élasticité, de résistance à l’humidité et ses capacités de 
flottaison, le polyuréthane est un matériau utilisé pour certaines planches (surf, planche à 

voile), parfois même sous la forme de sandwich peau compacte : cœur en mousse afin de 
donner de la flottabilité au matériel nautique ; 
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-tissus : utilisées sous forme de copolymères (par exemple polyuréthane : thermoplastique 
élastomère), les fibres sont drapées ou compactées pour la fabrication de vêtement ce qui leur 
confère de bonnes propriétés de confort (élasticité) et de respiration ; 

-industrie automobile : les mousses souples de polyuréthane sont employées dans les sièges, 
les appuie-tête, les accoudoirs, les mousses rigides de polyuréthane dans la volants, les toits 
ou les tableaux de bord ; 

-habitation/ chantiers : la mousse de polyuréthane associée a d’autres composants (panneaux 

de plâtre, de bois, de PVC.) permet de fabriquer une gamme de produits finis performants 

(panneau de toiture, panneau de doublage (isolation de mur), panneau de rénovation pour 
marché agricole. 

La fabrication de mousses polyuréthane conduit en général soit à un bloc, un boudin ou des 
feuillets de polyuréthane (lesquels doivent ensuit être découpés à la forme souhaitée), soit à 
une pièce directement moulée à la géométrie désirée. 

Dans le premier cas, le marché s’organise autour de : 

-la revente de blocs ou de boudins à des transformateurs ; 

- la transformation directe après moussage (ameublement, industrie automobile) ; 

- la revente de rouleaux (literie) ; 

- la découpe de plaques (literie) ; 

- la découpe numérique 3D (automobile). 

Dans le second cas, la mousse (ou le polyuréthane en cour de moussage) peut être injectée 
directement à l’intérieur de l’objet qu’elle est sensée garnir : on parle de moussage « in situ » ; 
par exemple : certains sièges automobiles sont fabriqués ainsi, le polyuréthane en cours de 
moussage est alors injecté directement dans la coiffe du siège, ce qui permet une meilleure 
adhésion mousse/ coiffe. Ce procédé est également adapté lorsque la géométrie de l’objet à 
garnir ne permet pas une découpe simple de blocs de mousse, comme dans le cas de 
l’isolation thermique dans les réfrigérateurs, congélateurs, etc. 
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POLVETHER 
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Fig.1.5 : Applications des mousses de polyuréthane flexible [42] . 


1.5.2. La chimie de polyuréthane. 


Les polyuréthanes résultent de la réaction de polymérisation chimique d’un isocyanate avec 
des groupements portant un hydrogène mobile (principalement des groupements hydroxyle), 
par exemple des fonctions alcool. Pour obtenir des matériaux alvéolaires, il faut coupler cette 
réaction à un dégagement gazeux qui permet la création d’alvéoles de gaz au sein du 
polymère [46] . 


Les polyuréthanes (PU) sont des polymères contenant des groupes uréthane dans la chaine 
moléculaire. Ils sont les éléments réactifs caractérisés par la présence du groupement 
(N=C=0). Les uréthanes sont produits par la réaction de polycondensation entre un 
isocyanate et un composé de polyol (hydroxyle). La réaction de polymérisation dans le 
l’uréthane peut être représentée par les équations suivantes 1511 : 


Une première réaction d’un isocyanate avec un alcool donne l’uréthane, cette réaction 
exothermique est connue sous le nom réaction de gélification. 
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R-N=C=0 + Ri —OH -►R-NH-C—0-R| 

II 

O 

Isocyanate Polyol Uréthane 


La réaction alcool-isocyanate peut être catalysée par une base de Lewis, en particulier les 
amines tertiaires. Les composés métalliques, et plus précisément les sels d’étain possèdent 
également un effet catalytique, cependant en tant qu’acide de Lewis, ils sont en général de 
plus faibles catalyseurs que les bases [52] . 


Les groupes uréthanes formés contient d’autre atomes d’hydrogène qui peuvent réagir ensuite 
avec un autre groupe d’isocyanate en transformation additionnelle pour donner des 
allophanates [53] . 


R—NH—C—O—R, + 
II 

O 


Uréthane 


R-N=C = Q 


Isocyanate 


O 

II 

R—N — C— O— CH - R, 

I 

R—NH—C 


Allophanate 


Une deuxième réaction de deux étapes dite réaction de soufflage ; il s’agit d’une réaction de 
l’eau avec des isocyanates : 

-La première étape concerne de la réaction entre l’isocyanate et l’eau, une réaction 
exothermique lente qui nécessite l’utilisation de catalyseurs et d’agitation [45] . Il se forme une 
amine avec dégagement de dioxyde de carbone qui joue un rôle important dans le soufflage 
du liquide. En outre, l’exotherme de la réaction joue également un rôle important dans 
l’expansion du produit. 


R—N=C=0 + H 2 0 -► R-NH 2 -I- C0 2 

Isocyanate Amine 

-Une deuxième étape est la réaction entre l’amine et un autre groupe isocyanate pour former 
une liaison urée di-substituée (figure 1.6), une réaction parasite entre l’urée di-substituée et un 
isocyanate peut se produire et conduire a la formation d’un biuret substitué figure (1.7). 
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R-N=C = 0 + RNH -► Ri -N-C-N-R 

H 6 H 

Isocyanate Amine Urée di-substituée 

Fig. 1.6 : Formation d’une Urée di-substitué. 


R-N=C = 0 
Isocyanate 


R,-N-C-N-R -► 

i h 

H O H 

Urce di-substituce 


O 

II 

R,-N-C-N-R 
C=0 

i 

N-R 

H 

bi u re t 


Fig.1.7 : Formation d’un biuret substitué. 

I.5.2.I. Autres réactions de l’isocyanate. 

Nous avons vu que les isocyanates réagissent préférentiellement avec les composés à 
hydrogène labile, mais d’autre réactions sont envisageables [54] . Les isocyanates peuvent 
réagir entre eux et conduire à la formation de dimères, trimères ou polymères en présence de 
catalyseur basique. Elle mène à la formation d’urétidione. 


O 

II 

^C s 

2 R —N = C = 0 -- R-N" N -R 

C 

II 

O 

urétidione 


Deux molécules d’isocyanate peuvent réagir et donner des carbodiimides [53] . 


: C=0 


Ri-N=C=N-R, +C0 2 

carbodiimide 
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Ensuite il peut se former un polycabodiimide [54] . 


R-N 

II 

n R — N=C = N — R -- ^ 

R 

polycabodiimide 


Par ailleurs, une fonction carbodiimide peut réagir à nouveau avec un isocyanate et former des 
urétonimines substituées [54] . 


R — N = C = N — R ♦ R —N=C = 0 


carlxxliimide 


R N-R 

N * 

N-C 
I I 
C-N 

* S 

O R 


co 2 


urétonimine 


La structure généralement formée par la trimérisation est celle d’un isocyanurate tri-substitué 

[ 53 , 54 ] 


3 R - N = C = O 


O 


R-N 

I 


N-R 

I 



N 

I 

R 



isocyanurate 


I.5.2.2. Isocyanates couramment utilisées dans la synthèse de polyuréthane. 

Les isocyanates sont des matières premières nécessaires pour la synthèse des polyuréthanes. 
Ils ont des structures aliphatiques, ou également des structures aromatiques. 


composants 


formule 


abréviation 

4,4-diisocyanate de 
diphénylméthane 

jT 

OCN 

-CH 2- 

jl^ 

NCO 

4,4-MDI 
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2,2-diisocyanate de 
diphénylméthane 

2,4-diisocyanate de 
diphénylméthane 

2,4-diisocyanate de toluène 
2,6-diisocyanate de toluène 

Tétraméthylxylènediisocyanate 



N = C = 0 



rsi = c= o 


2,2-MDI 


2,4-MDI 


2,4-TDI 

2,6- TDI 


TMXDI 


Isophoronediisocyanate 


Dicyclohexylméthane 4,4 
diisocyanate 

Hexaméthylène diisocyanate 


Isophorone diisocyanate 



OCN-(CH 2 )^-NCO 



IPDI 


H12 MDI 


HDI 


IPDI 


[C'y co cc~i 

Tab.1.1: Isocyanate couramment utilisées dans la synthèse de polyuréthane . 

I.5.2.3. Les polyols. 

Deux types de polyols sont couramment utilisés dans l’industrie, les polyéthers et les 
polyesters [45] . 
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I.5.2.3.I. Les polyols polyéthers. 

Ils représentent actuellement 80 à 90% des polyols utilisés sur le marché [45] . Les 
polyuréthanes préparés ont une excellente résistance à l’hydrolyse ; cependant, leurs 
propriétés mécaniques sont beaucoup plus faibles que les matériaux préparés à base de 
polyester [45] . 


R-O 



H 


Fig. 1.8: Polyol polyéther. 

I.5.2.3.2. Les polyols polyesters . 

Les produits obtenus à base de polyesters possèdent des caractéristiques plus élevées qu’avec 
des polyols polyéthers, mais une moindre résistance à hydrolyse. Leur prix est également 
supérieur à celui des polyols polyéthers [45] . 



Ri 



Fig. 1.9: Polyol polyester. 

I.5.2.3.3. Autres types de polyols. 

Deux autres types de polyols peuvent être utilisés dans la fabrication de polyuréthane, le 
polycarbonate polyols (figure 1.10) et le dimethyldantoine polyols (figure 1.11) [59] . 

HO-pR — O — C — O-fc R —OH 

II 

O 

Fig. 1.10 : Polycarbonate polyols. 
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ch 3 


ch 3 — c — C — O 


H —(OC 2 H 4 ) x -O-CH 2 CH 2 -N 


N-CH 2 CH 2 —O-(C2H4O V H 



II 


O 


Fig. 1.11 : Dimethylhydantoin polyols. 


I.5.2.4. Adjuvants 


I.5.2.4.I. Agents gonflants 


Pour la formulation des gaz d’expansion, deux méthode différentes sont usitées [45] : 


-Les agents gonflants par voie physique. 

-Les agents gonflants par voie chimique. 

Voie physique. 

Les agents gonflants par voie physique sont des solvants de bas point d’ébullition qui sont 

inertes dans les réactions chimiques et s’évaporent sous l’effet de l’exotherme de la réaction 

[ 60 ] 

Chlorofluorocarbones CFC. 

Avant 1990 les CFC (tableau 1.2) étaient utilisés comme un agent gonflant [45] . 

Le CFC-11, était le meilleur agent d’expansion physique, que ce soit pour les mousses 
souples ou bien rigides. Ces produits présentent l’avantage d’être non toxiques, bonne 
compatibilité avec les ingrédients de la réaction de moussage et une faible masse moléculaire, 
mais en revanche Les CFC sont responsables de la dégradation de l’ozone, ce qui conduit à 
leur interdiction en 1987 [60] . 


Abréviation 


formule Point d’ébullition °C ODP 


CFC-11 (R-11) CC1 3 F 

CFC-12 (R-12) CCLF 2 


-29.8 


23.8 


1 


1 


ODP : potentiel de destruction d'ozone 


Tab. 1.2 : Types de CFCs utilisés comme agent gonflant [60] . 
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Chlorofluorocarbures hydrogéné HCFCs. 


Un autre agent gonflant alternatif au CFC est le HCFC, comme indiqué dans le tableau. Ces 
agents gonflants ont toujours un problème d’appauvrissement de la couche d’Ozone et ont été 
progressivement éliminés en 2003 [60] . 


Abréviation 

formule 

Point d’ébullition °C 

ODP 

HCFC-21 

CHCLF 

8.9 

<0.05 

HCFC-22 

CHC10F 2 

-40.8 

0.055 

HCFC-31 

CH 2 C1F 

-9.1 

0.05 

HCFC-123 

CHCI 2 CF 3 

27.9 

0.02 

HCFC-124 

CHCIFCF 3 

-11 

0.02 

HCFC-132a 

CH 2 CICF 3 

6.1 

<0.05 

HCFC-132b 

CH 2 CICCIF 2 

46.8 

<0.05 

HCFC-141b 

CCI 2 FCH 3 

32.1 

0.11 

HCFC-142b 

CCIF 2 CH 3 

-9.8 

0.06 


Tab. 1.3 : Types de HCFCs utilisés comme agent gonflant [60] . 

Hydrofluorocarbures HFCs. 

Les HFCs ils ont été présentés comme une alternative aux HFCs, comme d’agents 
d’expansion. Le tableau (1.4) montre des agents gonflants HFC tels que le HFC-245fa et le 
HFC-365mfc. Toutefois, la forte contribution au réchauffement climatique et le coût élevé des 
HFC limitent également leur application . 

Ces composés tableau (1.4) possèdent l’avantage d’être sans effet sur la couche d’ozone. 


Abréviation 

formule 

Point d’ébullition °C 

ODP 

HFC-134a 

CF 3 CH 2 

-26.2 

0 

HFC-152a 

CHF 2 CH 2 

-24.7 

0 

HFC-245ca 

CHF 2 CHF 2 CH 2 F 

25.4 

0 

HFC-245fa 

CF 3 CH 2 CHF 2 

15.3 

0 

HFC-245eb 

CF 3 CHFCH 2 F 

22.7 

0 

HFC-365mfc 

CF3CH2CF2CH3 

40.2 

0 

HFC-365mffm 

CF 3 CH 2 CH 2 CF 3 

24 

0 

HFC-356 

CF 3 CH 2 CH 2 CH 3 

24.6 

0 

HFC-227ea 

CF 3 CHFCF 3 

30 

0 


Tab .1.4: Types de HCFCs utilisés comme agent gonflant [60] . 
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Perfliiorocarbures PFC 

Certains PFCs ont été proposés comme agents d’expansion [63] . Ils ne sont pas réputés 
dangereux pour la couche d’ozone, mais en raison de leur potentiel élevé en matière de 
réchauffement planétaire, l’agence de protection de l’environnement EPA a interdit leur 
utilisation [64] . 


Abréviation 

formule 

Point d’ébullition °C 

ODP 

PF-5040 

O 

Ph 

U 

-2 

0 

PF-5050 

c 5 h 12 

30 

0 

PF-5060 

c 6 h 14 

56 

0 

PF-5070 

C 7 H ,6 

80 

0 


Tab .1.5: Types de PFCs proposés comme agents d’expansion [60] . 

Hvdrofluorocarbures Ethers HFCEs. 

Certains HFCEs [65, 66 ], sont indiqués dans le tableau (1.6). Les HFC possèdent l’avantage 

d’une meilleure compatibilité avec les polyols durant le processus de moussage. 

Abréviation 

formule 

Point d’ébullition °C ODP 

HFCE-245 

CF 3 CH 2 -O-CF 2 H 

27 

0 

HFCE-356 

CF 3 CHFCF 2 -O-CH 3 

58 

0 

HFCE-347mcf 

CF 3 CF 2 CH 2 -O-CF 2 H 

46 

0 

HFCE-347pcf 

CF 2 HCF 2 -O-CH 2 CF 3 

56 

0 

HFE-254pc** 

CHF 2 CF 2 -O-CH 3 

37 

0 

HFE-236pc** 

CHF 2 CF 2 -O-CHF 2 

28 

0 

HFE-347mcc** 

CF 3 CF 2 CF 2 -O-CH 3 

34 

0 


Tab.1.6 : Types de HFCEs utilisés comme agent gonflant [Ashida 2007]. 

CO? liquide. 

Le dioxyde de carbone liquéfié apparaît comme agent gonflant en 1980 [60] . Cet agent 
d’expansion a un ODB nul, et un faible coût, raison pour laquelle son utilisation est de plus en 
plus populaire [62] . 
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C5-Hydrocarbures. 

Ces produits sont sans effet sur la couche d’ozone, mais ils ont le gros inconvénient d’être 
inflammables [45] . Le tableau (1.7) montre la liste des Cs-hydrocarbures usuels comme un 
agent gonflant. 


nom 

formule 

Point d’ébullition °C 

ODP 

Pentane 

CH 3 (CH 2 ) 3 CH3 

35 

0 

Isopentane 

C 2 H 5 CH(CH 3 ) 2 

30 

0 

Cyclopentane 

(CH 2 ) 5 

50 

0 


Tab .1.7: Types de Q-Hydrocrbures utilisés comme agent gonflant [60] . 

Voie chimique. 

La réaction de l’eau avec l’isocyanate produit le dioxyde de carbone et de la chaleur. La 
diffusion de dioxyde de carbone dans le polyol et l’exothermicité de la réaction provoquent 
l’expansion rapide du produit [45] . 



Fig.1.12: Mécanisme de la réaction isocyanate-eau [60] . 


D’autre agents gonflants peuvent être utilisés tels que, les composés organiques énolisables et 
l’acide borique [60] . 
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Typical 

Type Compounds Chemical Structure Enol Form 


C = N — OH Nitroalkane 
C = N —OH Aldoxime 


N = N —OH Nitrourea 
N = C — OH Acid amide 


Active methylene 
C = C—OH compounds 


R — CH 2 — NO, 

R — CH, — CH = N — OH 

H,N — CO — NH — NO, 

H R 

I I 

HN — C = 0 

R —CO—CH, —CO —R' 


O 

I 

R —CH = N —OH 
R — CH = CH — NH — OH 
O 

I 

H,N —CO —N = N —OH 
R 

I 

HN = C — OH 

R' 

I 

R — CO — CH = C — OH 


Tab.1.8 : Composées énolisables utilisés comme agent gonflant [60] . 



H 

I 



1 

—A = 

:B-OH 

Addition 


<j> —NCO 


\ Hs 
~ r 


Prototropy 

A = B - OCO — NH— <J> 

and 


v - 

décomposition 

A 


A = B + HO — CO — NH— <J> 

A 

C0 2 + H 2 N-«t> 

Décomposition 

| <j>—NCO 

<j) —NH —CO — NH —<{) 


Fig.1.13: Mécanisme de production de CO2 à partir des composées énolisables [60] . 
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A. Boronamidc formation 


HO-B-OH ♦ OCX'— 0 —NCO 
I 

OH 


—B—NH—0 — ♦ C0 2 | 
NH—♦ - 


B. Borate and Mater formation 


AAAA OH « HO-B-OH 
I 

OH 


\y\z\y\ O—B—o/vwx ♦ 3H 2 0 
I 

O 


O -NCO 


O-X’H-CO-NH-0 ♦ C0 2 


/\/\/\/\O O/VW\O-CO-NH-0-NH—CO-O /V/W-CO-NH-*- 


/ \ 


AAAAO 0/\/\/\/\0-CO-NH-0-N -CO -Q—/V\/X/\/X/X/VVX/\ 


I 


CO- NH- « - NH -CO —O —/w- 


Fig. 1.14: Mécanisme de la production de CO 2 à partir de l’acide borique [60] . 


I.5.2.4.2. Les catalyseurs [54 ’ 45] , 

La catalyse se réfère à l'accélération de la réaction, les acides tertiaires sont les catalyseurs 
basiques les plus utilisés. L’amine la plus couramment employée est 1,4- diazabicyclo (2,2,2)- 
octane (Dabco étant une des dénominations commerciales). 


On peut utiliser également les sels d’étain. L’octoate stanneux et le dibutyldilaurate d’étain 
(OBTDL) est le catalyseur organométallique le plus utilisé dans la fabrication de 
polyuréthane. 

I.5.2.4.3. Surfactants et agents de réticulation [45] . 

Les surfactants sont en général à base de polysiloxanes (silicones), ils jouent un rôle de 
stabilisateur physicochimique et de régulateur de cellules. 

Les agents de réticulation (allongeurs de chaîne) sont généralement des molécules à chaine 
courte contenant des amines ou fonction hydroxyle (les plus courant sont : le 1,4- butanediol 
et l’éthylène glycol, ou le glycérine). Ils peuvent être ajoutés pour fournir une augmentation 
de la charge ou la stabilité du produit final. 
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1.5.3. Méthodes de formulation des polyuréthanes cellulaires. 

La formation de polyuréthane cellulaire peut être réalisée par trois méthodes : la méthode 
« One shot », la méthode pré-polymère et la méthode quasi pré-polymère [46 ' 49 ' 50] . 

1.5 .3.1 Méthode de « One shot ». 

La méthode de «One shot » consiste à mélanger tous les composants (polyols, isocyanates, 
catalyseur, eau, solvant volatils) en même temps. Les inconvénients de cette méthode sont : le 
risque de dégradation du produit au cours de la réaction simultanée par l’effet de la 
température élevée, car cette réaction est exothermique [48] , en plus, le mélange de tous les 
composants au même temps peut s’avérer difficile. L’avantage de la méthode « One shot » est 
de préparer des grandes quantités avec un faible coût. 

1.5 .3.2. Méthode de pré-polymère. 

La méthode de pré-polymère consiste à faire réagir le polyol avec un excès d’isocyanate. 
Après la production de pré-polymère, les autres additifs, les catalyseurs et les agents 
moussants sont ajoutés. Cette méthode permet de contrôler facilement les propriétés 
physicochimiques du polymère fabriqué en changeant les quantités des additifs. 

1.5 .3.3. Méthode de Quasi pré-polymère. 

Comme pour la méthode pré-polymère la méthode Quasi pré-polymère consiste à mélanger le 
polyol avec l’isocyanate, à la différence qu’en utilisant au début seulement une partie de 
polyol ; Ensuite on provoque l’expansion par l’addition du polyol restant qui est déjà mélangé 
aux catalyseurs, aux agents gonflants et aux additifs. 

1.5.4. Méthodes de fabrication des polyuréthanes cellulaires. 

Deux voies sont possibles pour la mise en œuvre par moussage chimique, la coulée et 
l’injection [46] . 

I.5.4.I. Coulée souple. 

Ce procédé continu vise à fabriquer de grands blocs de mousse (mousse en bloc). Le mélange 
liquide, dont l’expansion commence très vite, est déposé en continu par une tète de coulée sur 
une bande transporteuse (figure 1.15). La bande de mousse peut facilement et rapidement 
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atteindre plusieurs mètres de largeur et de hauteur, tout en conservant une assez bonne 
homogénéité de la densité et de la taille des cellules. Tout au long de la coulée, la bande de 
mousse est entrainée a une vitesse de plusieurs mètres par minutes, avant d’être découpée en 
blocs. Les blocs sont ensuite stockés pendant environ 24 heures, temps nécessaire au 
polyuréthane pour achever sa polymérisation. Après la période de mûrissement, les blocs de 
mousse sont soit encore découpés en blocs plus petits soit refendus en feuilles et 
conditionnées en rouleaux. 



Fig. 1.15: Procédé mousse en bloc (Slabstock) [42] . 


I.5.4.2. Coulée par injection (expansion in-situ). 

A la différence de la coulée souple, ce procédé est séquentiel : ouverture de moule, injection, 
fermeture du moule, expansion, démoulage. Il vise à fabriquer de petites pièces de géométries 
bien spécifiques ou complexes. Une quantité très minutieusement déterminée est déposée par 
une tête d’injection dans le fond du moule (le trajet de dépôt peut même être optimisé pour 
améliorer la qualité de la mousse finale). Le mélange, dont l’expansion commence très vite, 
est confiné puisque le moule est fermé. Des évents doivent être placés sur le moule afin de 
permettre à l’air initialement présent qui y est de s’en échapper. 

L’objectif de cette opération est d’une part de remplir toutes les parties du moule, et d’autre 
part d’obtenir une homogénéité aussi parfaite que possible en terme de densité et de taille de 
cellules. Les mousses sont démoulées dés que le polyuréthane a une tenue mécanique 
suffisante, puis les blocs sont ensuite stockés pour le mûrissement (environ 24h également). 
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1.6. Conclusion. 

Ce premier chapitre a permis d’établir une bibliographie sur les polymères cellulaires, 
notamment le polyuréthane cellulaire, ses propriétés, ses fonctionnalités. Les bases de chimie 
de polyuréthane et les différentes méthodes de sa fabrication ont été présentées. 
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Analyse dimensionnelle d’une structure cellulaire 


2.1. Introduction. 


Le polyuréthane est un polymère compressible à cellules ouvertes qui est utilisé dans la 
fabrication des coussins pour fauteuil roulant pour des raisons de confort et de souplesse, du 
fait qu’il a une bonne capacité d’aération grâce à sa porosité, ainsi que sa capacité de la 
répartition du poids du corps. 

La mousse de polyuréthane peut subir des grandes déformations élastiques qui sont 
réversibles [36 ' 37] . 

La figure (2.1) : représente les réponses typiques en compression de trois types de mousses : 
(a) élastomère (caoutchouc, polymère, par exemple), (b) élastoplastique (métal) et (c) fragile 
(céramique) [21] . Lors d’une compression, le régime linéaire élastique est suivi pour tous ces 
matériaux d’un régime où la contrainte reste à peu prés constante, menant à un dernier stade 
où elle monte en flèche. Chacun de ces régimes correspond à un mécanisme de déformation 
typique. Au départ, les parois de la cellule subissent une flexion, ce qui donne un 
comportement élastique linéaire. Dès qu’une valeur critique est atteinte, les parois ou les 
arêtes commencent à s’effondrer. Dans le cas d’un matériau élastomère, les parois flambent 
de manière réversible. Pour les matériaux élastoplastiques, les parois se replient et, pour les 
matériaux fragiles, on observe la rupture des parois cellulaires aux endroits où les contraintes 
sont maximales. Finalement, les parois cellulaires opposées se touchent ou bien les parois ou 
arêtes commencent à s’entasser. C’est cela qui provoque la montée finale des contraintes en 
compression. 


8 

LU 

a. 

H 

LO 


(a) 



(b) 



(c) 


Fig.2.1 : Courbes contrainte-déformation en compression de trois types de mousses [34] . 
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2.2. Traction uniaxiale. 

Les courbes typiques contrainte-déformation des mousses en traction sont illustrées par la 
figure (2.2). D’après Gibson et Ashby [34] , la courbe de traction est assez différente du cas de 
la compression pour les trois types de matériaux ; seule la déformation initiale, qui est définie 
principalement par le fléchissement des parois et arêtes, en est similaire. Evidemment, le 
flambage ne joue aucun rôle pour un matériau cellulaire élastomère sous tension. Au 
contraire, les parois et arêtes de la cellule se tournent vers l’axe de tension et la rendent plus 
rigide. Les matériaux cellulaires plastiques se comportent de la même manière sous traction et 
sous compression, par la formation de zone plastiques, mais la rupture se produit pour des 
extensions beaucoup plus faibles en traction. La défaillance abrupte des matériaux cellulaires 
fragiles est observée à une valeur très inférieure à la résistance observée sous compression. 



Fig.2.2 : Courbes contrainte-déformation en traction de trois types de mousses [34] . 


On peut distinguer deux méthodes pour modéliser le comportement mécanique des solides 
cellulaires. Une méthode analytique, basée sur l’analyse dimensionnelle de structures 
périodiques et une méthode numérique, basée sur l’étude par éléments finis, en général elle se 
fait à l’aide d’un code de calcul qui traite ce genre de problème. 


2.3. Analyse dimensionnelle des solides cellulaires. 


Plusieurs travaux ont été effectués concernant l’analyse dimensionnelle des matériaux 
cellulaires : Hilyard [10] , ould Eyih [12] , Pradel [13] ; la majorité de ces études présente les 
mousses comme un assemblage de poutres représente les arêtes des cellules, avec des faces 
dans le cas des mousses à cellules fermées. 
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2.3.1. Nid d’abeilles. 

De manière générale, le paramètre le plus important qui contrôle les propriétés mécaniques est 
la densité relative [34] : 


O* 

d r = — =1-0 
P s 


( 2 . 1 ) 


Où : 


p*: est la masse volumique de la mousse. 

2 p s : est la masse volumique de matériau massif. 
O : la porosité de la mousse. 


Le nid d’abeilles est un arrangement régulier de cellules prismatique, comme le montre la 
figure (2.3). 


2-1 cos 0 —H 





Pour une charge appliquée suivant la direction xi. la théorie des poutres donne la flexion de la 
paroi : 


27 





Chapitre 2 


Analyse dimensionnelle d’une structure cellulaire 



Fig.2.4: charge appliquée suivant la direction xi. 


8 = Pr sin0/12E s I 


( 2 . 2 ) 


r 0 = 30°, dans un rayon d’abeille régulier. 

'i 

Où: I : le moment d’inertie de la paroi cellulaire (bt /12). 

^ E s : le module de Young du matériau cellulaire. 
b ; la profondeur de la paroi. 
v P : la charge appliquée suivant xi. 

La contrainte suivant la direction xi lié avec la charge appliquée P par la relation suivante : 


(Ti =- 

1 (h+lsinO) b 

La composante de la flexion donne la déformation dans la direction xi 

SsinQ 


(2.3) 


lcos9 


(2.4) 


Le module de Young parallèle à xi : 


E* 2 = 

£i 


(2.5) 


La combinaison de ces trois dernières équations donne le module de Young suivant la 
direction xi : 


HD 3 


cosd 


(h/l+sin9)sin 2 9 


( 2 . 6 ) 


Pour un hexagone régulier (h/1=1, 0= 30°) 
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Pour une charge appliquée suivant la direction x 2 : 


(2.7) 



t t t t t 

cr 2 


Fig.2.5 : charge appliquée suivant la direction x 2 . 


5=Wl 3 cos6/ 12E S I 

Où : (W: la charge appliquée suivant x 2 . 


La composante de la flexion donne la déformation dans la direction x 2 : 

Scos9 

£ 2 = -- 

h+l sin 9 

Le module de Young parallèle à x 2 : 



( 2 . 8 ) 


(2.9) 


( 2 . 10 ) 


La combinaison de ces trois dernières équations donne le module de Young suivant la 
direction^ : 


e 2 (t) 3 h/l+sin 9 
E s cos 3 9 

Pour un hexagone régulier (h/1=1, 0= 30°) 



( 2 . 11 ) 


( 2 . 12 ) 
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Figure.2.6 :Le module de Young théorique en fontion de module de Young experimental (a) 
nid d’abeille caoutchou (b) nid d’abeille métallique [38] . 


La densité relative d’un rayon d’abeille hexagonal : 

p* t h/l+2 

r PS l 2 cos6[h/l+sind) 


(2.13) 


Pour un hexagone régulier : 




2 

Vf 


t 

J. 


Ce qui donne (tenant compte de 2.7) : 



Les coefficients de Poisson pour des contraintes dans le plan est : 


(2.14) 


(2.15) 


e 2 cos 2 0 

e-l ( h/l+sind)sinO 


(2.16) 


v 


21 — 


£i 

e 2 


( h/l+sinO)sinO 
cos 2 9 


(2.17) 
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Pour un hexagone régulier : 


V 12~ V 21 - 1 


(2.18) 


Le module de cisaillement peut être calculé par la relation : 


G i 2 /t\ 3 ( h/f+sinO) 

E s ~ \{J ( h/£) 2 (l + 2h/l)cos8 


(2.19) 


Pour un hexagone régulier : 


E s 


m 3 


V3 Yt. 



( 2 . 20 ) 



2 


b 


Fig.2.7 : La variation du coefficient de poisson (a) et le module de cisaillement (b) pour un nid 

d’abeille de caoutchouc [38] . 


2.3.I.2. Comportement non élastique. 


L ’effondrement élastique 


Ce flambement apparait dans la phase plateau de la courbe contrainte-déformation, du a un 
flambement élastique des arêtes des cellules (figure 2.8). 
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ELASTfC 

BUCKLING 



Fig. 2.8 : Flambement élastique des arêtes. 


La force critique de flambage d’une poutre encastrée s’écrit : 


j-, _ n 2 n 2 E s I 

tcr ~ h 2 


( 2 . 21 ) 


Où : {n: facteur décrit le degré de la contrainte aux extrémités de la poutre. 
La contrainte élastique d’effondrement est donnée par la relation suivante : 


(°eÙ2- 


2blcos9 


( 2 . 22 ) 


La combinaison de ces deux dernières équations donne la relation suivante 


n 2 77-2 4-3 -I 

/ * \ _ n 7T t 1 rr 

°el 2~ 24 m 2 cosg 


(2.23) 


Pour un hexagone régulier : 


(a* el ) 2 =0.22(j) 3 E s 


(2.24) 


L ’effondrement plastique. 


Au passage de la force critique (le seuil de plasticité), les arrêtes flambent de manière 
plastique ; ce phénomène se traduit par des charnières aux jonctions des arrêtes, comme le 
montre la figure (2.9). 


32 



Chapitre 2 


Analyse dimensionnelle d’une structure cellulaire 


plastic 



Fig. 2.9 : L’apparition des rotules plastique. 


La contrainte d’effondrement plastique <j* t est donnée par les relations 


Suivant X 


■ . (v)i 

1 ' cr V ç Vt/ 


2 (h/l+sind) sinB 


Suivant X 2 : — l - 

cr v ? \i/ 2cos 


2 (9 


Pour un hexagone régulier : 

(v), _ (y) 2 _ 2 m 2 

s Gys 3 V l J 

Où : { o ys : la limite élastique du solide. 

La contrainte de cisaillement plastique s’écrit : 

(y) 12 _ ét\ 2 i 

Oys Vi/ (4 h/l)cos6 

Pour un hexagone régulier : 

(y) 12 _ _j_ fo 2 

(Tys 2 \ ; 3 \l J 

L ’effondrement par rupture. 

Dans cette étape on observe que les arrêtes se cassent figure (2.10). 


(2.25) 

(2.26) 


(2.27) 


(2.28) 


(2.29) 
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Dans cette partie, nous considérons h=l, toutes les équations développées dans ce cas sont 
pour des nids d’abeilles avec cette géométrie. 

Le module de Young suivant la direction X 3 est [14] : 


E l _ t 3 _ p* 

e s l 2cos0(l+sin0) ps 


(2.33) 


Les coefficients de Poisson v 21 = vf z sont supposés égaux au coefficient v s du matériau 
constituant le nid d’abeille. 


La relation de réciprocité permet d’obtenir les deux autres coefficients de Poisson : 


v 


13- 


E 1 

^ Vs ’ 


v- 


23- 


e 2 

t 3 Vs 


(2.34) 


Le comportement du nid d’abeille en cisaillement dépend de la géométrie. L’énergie de 
déformation a été exploitée pour estimer les modules de cisaillements G^et G 22 en fonction 
de la géométrie et la densité relative ainsi que G s le module de cisaillement du matériau 
constitutif [Zhang et Ashby 1992]. Les résultats sont : 


G 31 

G S 



(2.35) 


(l+sinÉ?) 2 f p*\ . Gî ? . 1 +sin 2 d ( p*\ 

4 (ps) ~ G s ~ 2 (ps) 


(2.36) 


Après la réponse élastique de la structure, un régime non linéaire sera suivi. Lorsqu’un nid 
d’abeille est comprimé suivant l’axe X 3 , et la contrainte dépasse le seuil d’élasticité, le 
flambage des parois se produit (figure 2 . 11 ) ; ce flambage caractérise la phase plateau dans la 
courbe contrainte-déformation. La charge critique de flambement est déterminée par la 
relation : 


K Es t 3 

F cr~ T\ T\ ~~T • 

(l-v|) l 

[K : facteur décrit la contrainte aux extrémités. 
Ou : < 

(/:la longueur des parois cellulaires. 

Lorsque b > 31, le constant K, est indépendant de b. 


(2.37) 
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b : la hauteur du nid d’abeille. 


La contrainte d’effondrement élastique est donnée par Zhang et Ashby 


[14] 


( °el ) 3 — 


10 F r - 


SK Ect* 


2cos6(l+sin6)l 2 [l-v 2 )cosd(l+sind)l 3 


(2.38) 


En substituant l’équation (2.33), v s = OA et K = 5.73 de Roark et Young [Roark et Young 
1975] pour satisfaire les conditions aux limites dans l’équation on obtiendra : 

{a* el ) 3 = 10.10 cos 2 0(l + sinG) 2 E s (yf (2.39) 

Noter que la fonction : cos 2 6( 1 + sinO ) 2 est presque constante si l’angle 0 est compris entre 
0° et 45°, tandis que la contrainte d’effondrement élastique est très sensible à la densité 
relative. 


Donc l’équation (2.39) se réduit à : 

<°A>3= 14-2 E s ($ 


(2.40) 


La contrainte d’effondrement plastique pour un hexagone régulier suivant X 3 est donnée par la 
relation [Wierzbicki 1983] : 


M 3 = 66 °ys ( 7 ) ' 

La contrainte à la rupture suivant X 3 s’écrit [Gibson et Ashby 1997]: 

/ * \ _ t _3_ p* 

[a cr>3- Ofsï 2cosg(1+sinQ) -°fs ps 


(2.41) 


(2.42) 


2.3.2. Modèle cubique [34] . 

Gent et Thomas [2 ' 3 ' 4] et Gent et Rusch r5 ' 6] furent les premiers modèles à réseau cubique, par 
la suite Gibson et Ashby [30] proposent un modèle cubique en trois dimension figure (2.12). 


Pour des mousses à faible densité (Z » t) l’analyse dimensionnelle, donne les relations 
suivantes : 


Mousse à cellules ouvertes : — oc 

Ps 



(2.43) 
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Mousses à cellules fermées : — oc 

Ps 


(2.44) 


Où : 


£ : Longueur des arêtes. 

J t : L’épaisseur des arêtes. 

t p : L’épaisseur des parois, dans le cas des mousses fermées. 



,c*tl tfly 



op*n <«ll facr I 




l 

ir 

i 


a 


Figure.2.12 : Cellule cubique proposée par Gibson et Ashby (a) ouverte (b) fermée [30] . 

2.3.2.I. Comportement élastique. 

Dans le régime élastique linéaire, Gibson et Ashby [30] considèrent que la mousse se déforme 
par la flexion des arrêtes horizontales figure (2.13). 
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Fig.2.13: La flexion des arrêtes horizontales dans le régime élastique linéaire [30] . 
Dans ce cas on a : 


F=oc a-8 2 


Où: fa : la contrainte de la flexion. 

F : la force exercée. 

La contrainte s’écrit : 


p 

G OC — 

C 


La déformation s’écrit : z ce - 


donc z oc — 

Eçl 


Sachant que : 5 oc 


pt 3 


Où : 1= — : le moment d’inertie d’une poutre. 


Le module de Young peut alors être calculé par : 


E* = - oc 


r £ n4 

1 


E r a 2 

ou — 

Eç 


(2.45) 


(2.46) 

(2.47) 

(2.48) 


(2.49) 


Pour le module de cisaillement G*, et le coefficient de poisson v* d’une mousse : 
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G* 


= C 2 d r 2 et:v*=e--l=C„. 


2C 7 


(2.50) 


Où : 


C 1 et C 2 sont des paramètres qui dépendent de la géométrie des cellules 


Ci = 1 et C 2 = - 

i / 8 


Nous obtenons alors : 


G* 3,2 * i 

- = -d r , v=- 

E s 8 ' 3 


(2.51) 


Dans le cas des mousses à cellules fermée, (figure 2.12.b) les relations sont les suivantes : 


E*= (p 2 E s d r 2 + (1- <p) E s d r + p( i 1 _^* ) 


(2.52) 


G*= - E s [(p 2 d 2 + (1 — (p)d r \ 
8 


, v =- 


(2.53) 


Où : 


p : La pression du gaz dans les cellules. 

(p : La fraction volumique de solide contenu dans les arêtes des cellules. 


Elle peut être calculée à partir de : 


Où : 


t 2 

^ t 2 +{Qlp)t p £ 

Q : est le nombre de faces par sommet, 
j P : le nombre moyen d’arêtes par faces dans une cellule. 


(2.54) 
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23 . 2 . 2 . Comportement non élastique. 

Après la phase élastique linéaire, la courbe contrainte déformation devient non linéaire. 
Gibson et Ashby on décrit deux types de déformations non linéaires 

L’effondrement élastique. 

Dans ce cas, les courbes contrainte-déformation en compression admettent une phase plateau, 
dû à un flambement élastique des arêtes des cellules, Comme le montre la figure (2.16.a). 



Figure 2.16 : Flambage des arêtes vertical (a) et l’apparition de rotules plastique au niveau de 

la jonction des arêtes (b). 


L’effondrement élastique est caractérisé par la contrainte élastique de plateau <j' el . Pour 
analyser le flambement, la formule d’Euler est appliquée aux poutres des cellules ; selon cette 
formule, la force critique de flambage d’une poutre encastrée est : 


Fcrit ~ 


n 2 n 2 E s I 

7 


(2.55) 


Où : [n : facteur décrit le degré de contrainte aux extrémités de la poutre. 


La contrainte élastique d’effondrement, est donnée par la relation suivante : 


o, 


el 


OC 


Ecrit 

i 2 


oc 


E£ 

e 4 


Ce qui donne : 


(2.56) 


— oc d r 2 : Tenant compte de (2.43) et (2.47). (2.57) 

Es 

Gibson [39] trouve que le facteur de proportionnalité = 0.05 
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Outre la contribution des arêtes, une forte dépendance entre l’effondrement dans une mousse 
à cellules fermées et la pression du gaz emprisonné dans les cellules se crée. La contrainte 
d’effondrement élastique s’écrit : 

a* el = 0.05 E s d r 2 + (p - Patm ) (2-58) 

Où : { P atm : la pression atmosphérique. 

L ’effondrement plastique. 

Les mousses rigides (polymères rigides, métaux), s’effondre d’une manière plastique. Ce 
phénomène, comme le flambage élastique, donne un plateau dans la courbe contrainte 
déformation à la différence qu’il n’est plus réversible. 

Pour les mousses à cellules ouvertes, cette ruine plastique se traduit par des rotules aux 
jonctions des arrêtes, Comme le montre la figure (2.16.b). 

La contrainte d’effondrement plastique a* t du matériau cellulaire, est donnée par la relation 
suivante : 


^ oc [d r ] 15 (2.59) 

Gy S 

Gibson [Gibson 1988] trouve que le facteur de proportionnalité = 0.3 

Dans le cas des cellules fermées, il faux rajouter la contribution des parois et celle du gaz 
emprisonné dans les cellules. Gibson et Ashby [30] obtiennent alors la relation suivante : 

— = 0.3 [q)d r ] is + (1- <p) d r + (2.60) 

GyS Gy S 

Où : { o ys : la limite élastique du solide. 
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ELASTIC COLLAPSE 
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RELATIVE DENSITY 


(a) 



(b) 


Fig.2.17.a : La contrainte de flambage élastique en fonction de la densité relative des 

mousses. 

b : Évolution de la contrainte plastique de la mousse en fonction de la densité 

relative [38] . 


L ’effondrement par rupture. 

Les mousse fragiles (les céramiques, polymères fragiles), s’effondrent par le cassement des 
arêtes des cellules. La contrainte à la rupture a^ r s’écrit : 

- Cellules ouvertes : 

a* r = 0.65[d r ] 1,5 (2.61) 

-Cellules fermées : 

— = 0.65[(pd r ] is + (1- cp) d r (2.62) 

a fs 

Où : {< a fs : le modulede rupturedu solide. 
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Figure.2.18 : Evolution de la contrainte de rupture de la mousse par le module de rupture du 

solide en fonction de la densité relative. 

2.3.3. Modèle tétraédrique. 

La structure tétraédrique est un assemblage de quatre poutres faisant toutes deux à deux un 
angle de 109°, cette description géométrique a été utilisée par Warren et Kraynik [35] . Ces 
derniers obtiennent la densité et le module de Young ainsi que le module de cisaillement et le 
coefficient de Poisson pour ce type de structure, et les résultats sont les suivants : 
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A 

2CÏÏ 2 ' 


(2.63) 



Z: la longeur des poutres. 
A: la section droite. 


E* 33V3 I , 2 

T s ~—xï dr 


(2.64) 


Où I : le moment d’inertie de la poutre. 
Section circulaire 


I/A 2 = 1/4tt, £ = 0.91cZ r 2 . (2.65) 

Section triangulaire : 

I/A 2 = 1/6V3, — = 1.1 d r 2 . (2.66) 

Es 

Section plateau (forme triangulaire dont les arêtes sont des arcs de cercle orientés vers 
Vextérieures) : 


I/A 2 


20V3-117T 
6( 2V3—7l) 2 


E* 

Es 


= 1.53 d r 


Le module de cisaillement peut être calculé par la relation : 

G„ _ lll 2 AI+12I 2 
E s 120V3i 4 / + 20V3 l 6 A 


Le coefficient de poisson est: 


5l 4 A 2 —3l 2 IA—36I 2 
101 4 A 2 +931 2 IA+36I 2 


(2.67) 


( 2 . 68 ) 


(2.69) 


2.3.4. Le tétrakaidécaèdre [8,9] . 


Ce polyèdre possède 14 faces, 6 carrées et 8 hexagonales, comme le montre la figure (2.20), 
les tétrakaidécaèdres sont des polyèdres semi réguliers qui peuvent paver l’espace avec des 
cellules de même forme dont la surface entre les cellules soit minimale 
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Le module de Young E{ 00 suivant les arêtes de cube est : 

r,* _ 6V2E5 I 

10 °- L 4( 1+1 2I/Al 2 ) 

{ / : le moment d'inertie des arêtes. 

A\ la surface de la section droide. 
i. la longeur des arêtes. 

Le coefficient de Poisson est de la forme : 


(2.70) 


* O c (Al 2 -121\ 

= 05 tend 

On note que pour les faibles densités relatives, le terme l 2 /E s A soit négligeable [8 9] 


(2.71) 


Et comme : 


1 _ 1 /12 1 2 i 4 \ 

E* 00 “ 6V2 \E S A E s l) 


(2.72) 
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D’où : 


j-,*_6-\/2 Es I 

^ 100 ~ 1 

Le module de cisaillement est donné par la relation : 

1 _ 2V2 1 2 V2 P /8 E s I+G t J 

G{ 2 E s A 6 E S I \5 E s I+G t J. 


(2.73) 


(2.74) 


Où : 


J: le moment dipolaire des arêtes. 

[Gp.le module de cissaiment du matériau massif. 


2.4. L’énergie absorbée pour une mousse élastique. 

2.4.1. Construction des diagrammes contrainte-déformation [40] . 

La densité relative p(e)/ p s après une déformation compressive £ s’écrit : 

p( g) = £_(. 1 3 

Ps Ps Vl-e/ 


(2.75) 


Où : — est la densité relative initiale. 

Ps 

Lorsque la densification est complète, et la courbe contrainte-déformation devient presque 
verticale, où =1, la déformation devient : tf = 1 -p/p s . 

Maiti, Gibson et Ashby [40] trouvent expérimentalement que la fin du plateau correspond à : 
= - , (le solide occupe - du volume total), et la déformation devient : e s = 3 -p/p s . 

Ps 3 3 


La longueur des arêtes, après une certaine déformation £ peut s’écrire comme suit: 

1(0) - 1 \ _ V Ps ’ _[ 

U ; l ° 1 -(P/Ps ) l/3 


(2.76) 


Où : {l 0 : la longueur initiale des arêtes. 


En utilisant l’équation (2.75) on obtiendra : 


Ke) =l 0 


U-( P 1 ) l/ 3 l 

I Um > 

1 -(P/Ps ) l/3 • 


(2.77) 


En remplaçant F par la valeur a Ïq et tenant compte (2.21) en obtiendra 
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De même on obtiendra pour les mousses plastiques : 



(2.78) 


(2.79) 



Fig.2.21 : diagrammes contrainte-déformation [40] . 

2.4.2. Le phénomène de densification. 

L’application très fréquente des mousses polymériques est l’emballage de matériel fragile. 
Cette application est due à la capacité des mousses d’absorber de l’énergie en compression. 
Pour une description simplifiée, considérons une courbe typique contrainte-déformation d’un 
matériau cellulaire soumis à un essai de compression comme celle de la figure (2.22) : la 
courbe se décompose en trois zones : élasticité linéaire, plateau d’effondrement et 
densification. 

W : énergie absorbée jusqu'à la densification, 
w 

— : Air sous la courbe contrainte-déformation. 

E s 
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Fig 2.22: Courbe contrainte-déformation pour un matériau cellulaire. 

2.4. 3. Construction des diagrammes de l’énergie absorbée 

L’énergie absorbée jusqu'à une déformation £*, où la contrainte est a est donnée par : 

W= / £ °crd£ + f E aie)de 

U £q 



--eg + f d£ 

2 E s u J £o E s 


r£*a(e) 


(2.80) 

(2.81) 


On intègre cette équation en utilisant l’équation (2.78) pour les solides cellulaires élastiques, 
et l’équation (2.79) pour les solides cellulaires plastiques : on obtiendra les diagrammes : 




Fig.2.23 : les diagrammes de l’énergie absorbée. 
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2.4. 4. La courbe enveloppe. 

L’enveloppe de l’énergie absorbée est l’ensemble des points critiques associé à la famille des 
courbes (figure 2.23), c’est une courbe tangente à chacun des diagrammes pour différentes 
densités, elle permet d’estimer, en chaque section de densité, la valeur maximale de l’énergie 
qui peut être absorbée par le polymère cellulaire pour une contrainte bien déterminée. 

Pour les mousses élastiques l’équation de la droite correspondante à l’enveloppe est 
approximativement : 

l=o.n(£) /a (2.82) 

<7 V 

Pour les mousses plastiques avec — = 0.1, l’équation de la droite est approximativement : 

Es 

W /<x\ 2 / 3 

- = 0.05(-) (2.83) 

2.5. Conclusion. 

En conclusion, l’analyse dimensionnelle permet de trouver des expressions reliant les 
caractéristiques macroscopiques d’un matériau à sa micro structure. Nous avons pu voir que 
les caractérisations mécaniques d’une mousse sont fortement liées de sa structure et des 
propriétés du matériau constitutif. La plupart des modèles, considèrent les solides cellulaires 
comme des assemblages de polyèdres réguliers, forme cubique, tétraèdre, et tétrakaidécaèdre. 
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3.1. Introduction. 


La simulation numérique est l'un des outils permettant de simuler des phénomènes réels. 
Appelée aussi simulation informatique, elle désigne un procédé selon lequel on exécute un 
programme informatique sur un ordinateur en vue de simuler par exemple un phénomène 
physique complexe (ex: chute d'un corps sur un support mou, résistance d'une plateforme 
pétrolière à la houle, fatigue d'un matériau sous sollicitation vibratoire, usure d'un roulement à 
billes, ...). 

Les simulations numériques scientifiques reposent sur la mise en œuvre de modèles 
théoriques utilisant souvent la technique dite des éléments finis. Elles sont donc une 
adaptation aux moyens numériques de la modélisation mathématique, et servent à étudier le 
fonctionnement et les propriétés d'un système modélisé ainsi qu'à en prédire son évolution. On 
parle également de calcul numérique. 

Les interfaces graphiques permettent la visualisation des résultats des calculs par des images 
de synthèse. 

Dans la plupart des projets [23] , la MEF est souvent utilisée pour prédire la déformation de la 
structure en service et fournir des données relatives à la phase de construction (déplacement 
et déformation en cours de construction). Le calcul par la MEF est souvent effectué afin de 
vérifier le comportement constaté du dispositif. 

L'analyse de type éléments finis d'une structure peut notamment comporter les étapes 

[311 

suivantes : 

- la définition de la géométrie et conditions aux limites, 

- le choix de la relation déformations-contraintes des matériaux utilisés, 

- les caractéristiques mécaniques des interfaces, 

- la définition des chargements, 

- la visualisation des résultats. 
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3.2. Le code de calcul Abaqus 

3.2.1. Structure générale d’un code de calcul 

Pré Processeur 
Choisir le type d’élément 
Entrer les propriétés géométriques 
Entrer les propriétés physiques 
Créer le modèle géométrique 
Créer le maillage 
Appliquer les sollicitations 
Imposer les conditions aux limites 


_V_ 

Solveur 

Choisir le type d’analyse (dynamique, statique ...) 
Construire la matrice et les vecteurs élémentaires [Ke], [Fe] 
Assembler [Ke]et [Fe] dans [K], [F] 

Prendre en compte les conditions aux limites 
Résoudre le système d’équation [K]{U}=[F] 
Calculer les variations additionnelles 


V 

Post-Processeur 

Présenter les résultats de façon intelligible et synthétique : 

* Sous forme numérique 

* Sous forme graphique 

Effectuer des fonctions complémentaires : combinaisons, 
animation,... 

Fig.3.1 : La structure générale d’un code de calcul. 

3.2.2. Présentation générale d’Abaqus. 

ABAQUS est un code de calcul utilisant la MEF créé en 1978. Selon la méthode de 
résolution, on distingue deux grands codes [31] : 
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- ABAQUS/Standard : résolution par un algorithme statique implicite, 

-ABAQUS/Explicit : résolution par un algorithme dynamique explicite. 

La version que nous avons utilisé, ABAQUS 6.9.1, permet de traiter des problèmes 
couplés (hydromécanique, thermomécanique ...) dans les divers domaines (mécaniques, 
thermiques, électriques ...). Un très large choix d’éléments en 2D ou 3D est proposé dans le 
code pour des analyses statiques, dynamiques, linéaires et non-linéaires. La structure d'une 
analyse ABAQUS est présentée sur la figure (3.2). 


Pré processeur 



Fichier .inp, .cae 



Post-processeur 


Fig.3.2 : La structure d’une analyse ABAQUS [22] . 

Le pré-processeur lit le fichier d'entrée qui définit la géométrie du modèle, sa 
discrétisation, les conditions aux limites, les chargements appliqués et les résultats demandés. 
Ensuite ABAQUS exécute le calcul, le post-processeur se charge de récupérer les résultats et 
permet aux utilisateurs de les visualiser puis de les exploiter. Plus de détails concernant 
l'utilisation d'ABAQUS seront donnés par la suite [22] . 

- Domaines physiques: 

Mécanique, Thermique, Électrique, Problèmes couplés. 
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-Problèmes: Statique et dynamique, Linéaires et non linéaires. 

-Catégories d’éléments: Milieu continu (2D et 3D), Poutres, plaques, coques, éléments 
spéciaux (ressorts, masses,...) 

-Géométries: Approximations linéaires et quadratiques, triangles, quadrangles, hexaèdres, 

tétraèdres, pentaèdres 

-Un très large choix d’éléments. 

-Possibilités de programmer de nouveaux éléments (en FORTRAN, sur ABAQUS/Standard). 


3.2.3. Le modèle Hyperfoam Ogden sur Abaqus 

Le modèle « Hyperfoam » est destiné à la mousse élastique, ce modèle a été utilisé par 
plusieurs auteurs, Mills et al [24] , Mills et Lyn [25] , Ont utilisé le modèle Hyperfoam Ogden 
pour simuler le comportement mécanique des mousses de différentes applications, des 
semelles pour souliers et casque anti-choc. Grujicic et al [26] ont étudié la résistance en 
compression pour certains polymères cellulaires utilisés dans le secteur automobile et 
transports. Leur étude est basée sur le modèle Hyperfoam Ogden. 

3.2.3.I. L’équation constitutive générale pour les matériaux hyperélastique. 

Pour décrire le comportement viscoélastique à l’état compressible des matériaux, une 
équation constitutive viscoélastique est utilisée. Le tenseur de la contrainte S est supposée 
divisé en deux composantes. [27] : l’une d’équilibre, l’autre de surcharge. 

s = S G +S GV (3.1) 


La fonction énergie de déformation est donnée par formule : 

W = J S **D avec J= det F 


Où : 


w: la fonction energie de déformation. 
F\ le gradient de déformation. 

S\ le tenseur de contrainte de CAUCHY. 
D \ le tenseur de déformation. 


(3.2) 


On peut déterminer la composante du tenseur de déformation par la relation suivante : 



(3.3) 
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Avec : 


W= W(U) = W(C) 


> 0 pour C ^ I 
= 0 pour C = / 


(C\ le tenseur de la droite de CAUCHY — GREEN. 
[U‘. le tenseur des allongements. 


(3.4) 


Il est à noter que : (C=I) dans le cas non déformé, (C=É I) dans le cas déformé. 


L’équation constitutive générale pour le cas non-linéaire hyperélastique est obtenue en 
substituant les équations (3.3) et (3.4) dans (3.2) on aura : 


S = 2 J- 1 F- ?!m -F T 

dC 

3.2.3.2. La fonction énergie de déformation pour les polymères compressibles. 


(3.5) 


Pour décrire le comportement mécanique des polymères dans le cas compressible, la fonction 
énergie de déformation a été proposée par Hill [28] et Storakers [29] : 


N 


W 


-I 

i=i 


Ei_ 

a} 


î a i j_ ii a i i l° 
/ L-i I An I An 


3 + 


Pi 


f) 




- 1 ) 


(3.6) 


En substituant l’équation (3.6) dans l’équation (3.5) on obtient : 




k= 1 i= 1 


s = 2r LLia3 À ‘ ,+ tt l 1 ~âr nknk 


(3.7) 


(N\ l'ordre de calibration. 

Où : Pi , ai, et (3p. sont des constantes de calibration depond de la tepérature. 


Il est à noter que : / (/) = j J el ) — 1^ 

(3.8) 

x i = (j th )-lx i _ x 1 x 2 k=J el 

(3.9) 


{ J el \ est le rapport de volume lors d'une déformation élastique. 
J th \ est le rapport de volume lors d'une dilatation thermique. 
Xi'.sont les allongements principaux. 

J el : Est lié à ] th par l’équation (3.10) 
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J el =y- (3.io) 

J th 

Sachant que : J est le rapport de volume total. 

Le coefficient de Lamé initial est : 


N 

I 

i= 

Et le module de cisaillement initial est : 

N 

i =1 

Le coefficient // détermine le degré de compressibilité, il est relié au coefficient de poisson v L 
par l’expression suivante : 

Pi = i _ v ^ respectivement v t = 1= 1,2,., N (3.14) 

Pour le cas particulier : = /? = const, v : est le coefficient de Poisson. 

Selon Hill |2X| et Storakers [Storaker 1986] les considérations suivantes sont toujours valides. 

PiCCi >0 (i = 1,2,. ,N) et /? > —^ (3.14) 

3.3. Géométrie de la structure considérée et conditions aux limites . 

Il s’agit d’une structure cubique de dimensions 40*40*40mm, soumise à une compression 
uniaxiale et uniforme, Les déplacements et les rotations sont fixés sur le coté inferieur de la 
structure. La figure (3.3) résume l’ensemble de la modélisation. 


Pi = Po 


+ 1/3) = K 0 


(3.11) 


(3.12) 
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Fig. 3.3 : l’ensemble de la modélisation. 

3.4. Matériau considéré. 


Notre choix s’est porté sur les mousses polyuréthane à cellules ouvertes, Ces mousses 
couvrent des plages de densité apparente très larges, présentent des morphologies et structures 
très différentes, sont proposés avec des modules allant d’une grande souplesse jusqu'à la 
rigidité et couvrent la majorité des applications connues [44] . 

Le polyuréthane massif possède les caractéristiques suivantes [34] : 

Le module d’Young : Es = 2500 MPA. 

La densité : d s = 1.18g /cm 3 . 

La limite élastique : a ys = 40MPa 

3.4.1. Calcul des constantes élastiques. 

Le matériau cellulaire considéré est le polyuréthane a cellule ouverte forme cubique, les 
formules (2.49, 2.50, 2.51) développées par Gibson et Ashby, [30] , ont permis d’obten ir le 
module de Young E*, le module de cisaillement G* ainsi que le coefficient de poisson v*. 
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La densité relative d r 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

la densité d (g /cm 3 ) 

0.354 

0.472 

0.59 

0.708 

0.826 

0.944 

Module de Young E* 

(MPa) 

225 

400 

625 

900 

1225 

1600 

Module de 

cisaillement G*(MPa) 

84.4 

150 

234.4 

337.5 

459.4 

600 

Le coefficient de 

0.33 

0.33 

0.33 

0.33 

0.33 

0.33 


poisson v* 


Tab.3.1 : Calcul des constantes élastiques. 

3.5. Maillage de la structure. 

Le type d’élément adopté pour le maillage de la structure est l’élément tétraédrique à quatre 
nœuds. Ce choix est justifié par le fait que son application large dans les problèmes d’analyse 
de contrainte a conduit à des résultats très appréciables. Dans la figure 3.4 on montre un 
maillage composé de 39963 éléments. 



Fig. 3.4 : Maillage de la structure. 
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3,6. Résultats et discussion. 

3.6.1. Visualisation des résultats pour une densité 30%. 

Contrainte de Von Mises. 

Un système de forces extérieures appliqué à un corps engendre un système de forces 
intérieures résultant de l’interaction des particules élémentaires. A chaque force élémentaire 
correspond une contrainte. La contrainte de Von Mises est utilisée pour illustrer la répartition 
des niveaux de contrainte subis par le matériau. La contrainte de Von Mises s’écrit : 

— a yy) (°yy — °zz) + (°zz — °x ) 2 % 6 (t \y + r | 2 + Ty Z ) 

°vm = I-~- S <*e 

\ L 


Où les valeurs de a sont les contraintes normales et t les contraintes de cisaillement. 



S, Mises 
(Avg: 75<>/o) 

+ 2.105e + 01 
+ 1.967e+ 01 
+ 1.829e+ 01 
+ 1.691e+01 
+ 1.552e+ 01 
+ 1.414e+01 
+ 1.276e+ 01 
+ 1.138e+01 
+ 9.994e+O0 
+ 8.612e + O0 
+ 7.229e + O0 
+ 5.847e + O0 
+ 4.465e + O0 


Fig.3.5 : Contraintes de Von Mises pour d r = 30%. 


La déformation. 


Le polyuréthane cellulaire est un solide déformable. Au cours de sa compression. Elle se 
déforme. La grandeur des déformations est définit comme des élongations par unité de 
longueur dans une direction. La déformation s suivant une longueur L s’écrit : 


£ 


Al 

L 
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E, Max. Principal 
(Avg: 75°/o) 

+ 8.074e-01 
+ 7.401e-01 
+ 6.728e-01 
+ 6.055e-01 
+ 5.383e-01 
+ 4.710e-01 
+ 4.037e-01 
+ 3.364e-01 
+ 2.691e-01 
+ 2.019e-01 
+ 1.346e-01 
+ 6.731e-02 
+ 3.544e-05 


Fig.3.6 : Déformation pour d r = 30%. 


Courbe contrainte déformation. 


30 % 



Fig.3.7 : Réponse mécanique pour d r = 30%. 
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La figure (3.7) présente la réponse mécanique lors d’un essai de compression uniaxiale pour 
une valeur de densité relative 30%, selon cette courbe on peut distinguer trois zones. Dans la 
première zone, la contrainte augmente considérablement de manière quasi proportionnelle à la 
déformation, cette zone correspond à la réponse élastique. La deuxième zone est caractérisée 
par un net ralentissement de l’accroissement de la contrainte à partir d’environ 5% de 
déformation. La troisième zone se traduit par une nouvelle augmentation rapide de la 
contrainte qui caractérise le processus de la densification de la mousse. 


La densité 

Limite 

Déformation à 

Déformation 

Contrainte 

relative d r 

élastique a* el 

la 

maximale £ max 

maximale o max 

(30%) 

(MPa) 

densification e d 


(MPA) 

Valeur simulée 

11.38 

0.6 

0.8 

21.05 


Gibson et 11.25 

Ashby 

Ecart (%) 1.14 


Tab.3.2 : Résultats obtenus pour d r = 30%. 
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3.6.2. Comparaison entre l’énergie absorbée par le solide cellulaire et le solide qui la 

compose. 


Densité 50%. 



S, Mises 
(Avg: 75®/o) 

+ 3.919e + 01 
+ 3.650e + 01 
+ 3.382e + 01 
+ 3.114e+ 01 
+ 2.845e + 01 
+ 2.577e + 01 
+ 2.309e + 01 
+ 2.040e + 01 
+ 1.772e + 01 
+ 1.504e + 01 
+ 1.235e + 01 
+ 9.669e + 00 
+ 6.986e + 00 


Fig.3.8 : Contraintes de Von Mises pour d r = 50%. 



E, Max. Principal 
(Avg: 75o/o) 

+6.090e-01 
+ 5.582e-01 

- +5.075e-01 

- +4.567e-01 

- +4.060e-01 

- + 3.553e-01 

- +3.045e-01 

- + 2.538e-01 

- +2.030e-01 

- + 1.523e-01 

- + 1.015e-01 

- +5.077e-02 

- + 2.673e-05 


Fig.3.9 : Déformation pour d r = 50%. 
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Fig.3.10 : Comparaison entre la réponse mécanique du polyuréthane de d r = 50% et le 

polyuréthane massif. 


La figure (3.10) est un schéma où est représentée la réponse en compression d’un 
polyuréthane dense et d’une mousse fabriquée à partir de ce polyuréthane. 

On remarque que pour une même contrainte o* el - 31.14 MPa, le polyuréthane cellulaire 
absorbe une quantité d’énergie (les aires hachurées) beaucoup plus importante que celle 
absorbée par le solide qui la compose (le polyuréthane massif).ce qui implique que, pour la 
même quantité d’énergie absorbée, la mousse produit un moindre pic de contrainte. 

Comme la figure (3.10) le montre, c’est le long plateau de la courbe contrainte- déformation 
qui permet une grande absorbation d’énergie à contrainte a* t presque constante jusqu'à une 
déformation plus de 60%. 
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La densité 

Limite 

Déformation à 

Déformation 

Contrainte 

relative d r 

élastique a * L 

la 

maximale £ max 

maximale o max 

(50%) 

(MPa) 

densification e d 


(MPA) 

Valeur simulée 

31.14 

0.4 

0.6 

39.19 


Gibson et 

31.25 

Ashby 


Ecart (%) 

0.35 


Tab.3.3 : Résultats obtenus pour d r = 50%. 


3.6.3. Influence de la densité. 

Pour mettre en évidence l’influence de la densité sur la réponse mécanique des solides 
cellulaires de polyuréthane, la méthode de modélisation a été appliquée à des mousses de 
densité relative 40%, 70% et 80% 

Densité 40%. 



S, Mises 
(Avg: 75®/o) 

+ 3.184e + 01 
+ 2.966e + 01 
+ 2.748e+ 01 
+ 2.530e + 01 
+ 2.312e + 01 
+ 2.094e + 01 
+ 1.876e + 01 
+ 1.658e + 01 
+ 1.440e+ 01 
+ 1.222e + 01 
+ i.004e+ 01 
+ 7.856e + 00 
+ 5.676e + 00 


Fig.3.11 : Contraintes de Von Mises pour d r = 40%. 
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E, Max. Principal 
(Avg: 75®/4>) 

+ 6.920e-01 
+ 6.344e-01 
+ 5.767e-01 
+ 5.190e-01 
+ 4.614e-01 
+ 4.037e-01 
+ 3.460e-01 
+ 2.884e-01 
+ 2.307e-01 
+ 1.730e-01 
+ 1.154e-01 
+ 5.770e-02 
+ 3.038e-05 


Fig.3.12 : Déformation pour d r = 40%. 


Densité 70% 



S, Mises 
(Avg: 75o/o) 

+ 8.082e + 01 
+ 7.529e + 01 
+ 6.975e + 01 
+ 6.422e+01 
+ 5.868e + 01 
+ 5.315e + 01 
+ 4.761e+01 
+ 4.208e+ 01 
+ 3.655e + 01 
+ 3.101e + 01 
+ 2.548e + 01 
+ 1.994e + 01 
+ 1.441e + 01 




Fig.3.13 : Contraintes de Von Mises pour d r = 70%. 
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E, Max. Principal 
(Avg: 75<>/o) 

+ 4.308e-01 
+ 3.949e-01 
+ 3.590e-01 
+ 3.231e-01 
+ 2.872e-01 
+ 2.513e-01 
+ 2.154e-01 
+ 1.795e-01 
+ 1.436e-01 
+ 1.077e-01 
+ 7.181e-02 
+ 3.591e-02 
+ 1.891e-05 


Fig.3.14 : Déformation pour d r = 70%. 


Densité 80%. 



S, Mises 
(Avg: 75<>/o) 

+ 9.062e + 01 
+ 8.441e + 01 
+ 7.821e + 01 
+ 7.200e + 01 
+ 6.580e + 01 
+ 5.959e+ 01 
+ 5.339e+ 01 
+ 4.718e + 01 
+ 4.098e + 01 
+ 3.477e + 01 
+ 2.857e + 01 
+ 2.236e + 01 
+ 1.615e + 01 




Fig.3.15 : Contraintes de Von Mises pour d r = 80%. 


66 
































Chapitre 3 


simulation numérique 



E, Max. Principal 
(Avg: 75«>/o) 

_ +4.055e-01 
™ +3.717e-01 

- +3.379e-01 

- +3.041e-01 

- + 2.703e-01 

- + 2.365e-01 

- +2.028e-01 

- + 1.690e-01 

- + 1.352e-01 

- + i.0i4e-0i 

- +6.760e-02 

- +3.381e-02 

- + 1.780e-05 


Fig.3.16 : Déformation pour d r = 80%. 



Fig.3.17 : Réponse mécanique de polyuréthane de 3 densités relatives 40%, 70%, 80%. 
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Le comportement mécanique des mousses est étudié en fonction de leur densité relative. La 
figure (3.17) présente leur réponse mécanique lors d’un essai de compression uniaxiale. Les 
courbes sont tracées pour trois valeurs de la densité relative 40%, 70%, 80% sur le même 
graphique. Pour les mousses de densité relative 70% et 80%, on observe un pic de contrainte 
avant d’atteindre le plateau, ce qui signifie que le début de la densification se produit d’autant 
plus tôt que la densité est élevée. En effet, plus la mousse est dense, plus les vides sont 
comblés rapidement. 


La densité 

Limite 

Déformation à 

Déformation 

Contrainte 

relative d r 

élastique o* el 

la 

maximale £ max 

maximale o max 

(40%) 

(MPa) 

densification e d 


(MPA) 

Valeur simulée 

20.94 

0.52 

0.69 

31.84 


Gibson et 20 

Ashby 

Ecart (%) 4.49 


Tab.3.4 : Résultats obtenus pour d r = 40%. 


La densité 

Limite 

Déformation 

Contrainte 

relative d r 

élastique o* el 

maximale £ max 

maximale o max 

(70%) 

(MPa) 


(MPA) 

Valeur simulée 

/ 

0.43 

80.82 


Gibson et 61,25 

Ashby 


Tab.3.5 : Résultats obtenus pour d r = 70%. 
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La densité 

Limite 

Déformation 

Contrainte 

relative d r 

élastique a* el 

maximale £ max 

maximale o ma x 

(70%) 

(MPa) 


(MPA) 

Valeur simulée 

/ 

0.4 

90.62 


Gibson et 80 

Ashby 


Tab.3.6 : Résultats obtenus pour d r = 80%. 

3.6.4. Influence de trou. 

Les paramètres du polyuréthane cellulaires troué ont été les mêmes que ceux du solides 
cellulaire non troué ; pour pouvoir comparer les résultats, on choisit dans ce cas une densité 
de 60%. 


Maillage de la structure trouée. 




La figure (3.18) montre un maillage assez raffiné prés du trou composé de 46472 éléments 
dans le cas du trou horizontal contre 46206 éléments dans le cas vertical. 
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S, Mises 
(Avg: 75%) 

+ 6.123e+ 01 
+ 5.703e + 01 
+ 5.284e+ 01 
+ 4.865e+ 01 
+ 4.446e+ 01 
+ 4.026e+ 01 
+ 3.607e+ 01 
+ 3.188e+ 01 
+ 2.769e+ 01 
+ 2.349e+ 01 
+ l,930e+ 01 
+ 1.511e+01 
+ 1.092e+ 01 


Fig.3.19 : Contraintes de Von Mises pour d r = 60% (non trouée). 



E, Max. Principal 
(Avg: 75%) 

+ 5.046e-01 
+ 4.626e-01 
+ 4.205e-01 
+ 3.785e-01 
+ 3.364e-01 
+ 2.944e-01 
+ 2.523e-01 
+ 2.103e-01 
+ 1.682e-01 
+ 1.262e-01 
+ 8.412e-02 
+ 4.207e-02 
+ 2.215e-05 


Fig.3.20 : Déformation pour d r = 60% (non trouée). 
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«— 60 % 



Fig.3.21: Réponse mécanique pour d r = 60% (structure non trouée). 


La densité 

Limite 

Déformation à 

Déformation 

Contrainte 

relative d r 

élastique a * L 

la 

maximale £ max 

maximale o max 

(40%) 

(MPa) 

densification c D 


(MPA) 

Valeur simulée 

44.46 

0.3 

0.5 

61.23 


Gibson et 45 

Ashby 

Ecart (%) 1.2 


Tab.3.7: Résultats obtenus pour d r = 60%. (Structure non trouée). 
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ODB: Joo-7.odo Aoaqas/'Sca'rta'd 6. lZ-i Tnu 


S, Mises 
(Avg: 75<>/o) 

+ 6.014e+01 
+ 5.598e+ 01 
1 5.103e 1 01 
+ 4.767e+ 01 
+ 4.351e+ 01 
+ 3.935e+ 01 
+ 3.520e+ 01 
+ 3.104e+ 01 
+ 2.688e+ 01 
+ 2.273e+0i 
+ 1.857e+ 01 
+ 1.441e+ 01 
+ 1.025e+ 01 


Fig.3.22 : Contraintes de Von Mises pour d r = 60% (trou vertical). 



E, Max. Principal 
(Avg: 750/b) 

+ 5.058e-0i 
+ 4.636e-01 
+ 4.215e-01 
+ 3.793e-01 
+ 3.372e-01 
+ 2.950e-01 
+ 2.529e-01 
+ 2.107e-01 
+ 1.686e-01 
+ 1.264e-01 
+ 8.430e-02 
+ 4.215e-02 
+ 1.304e-06 


Fig.3.23 : Déformation pour d r = 60% (trou vertical). 
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S, Mises 
(Avg: 75o/o) 

+ 8.942e + 01 
+ 8.216e + 01 
+ 7.491e + 01 
+ 6.766e + 01 
+ 6.040e + 01 
+ 5.315e + 01 
+ 4.590e + 01 
+ 3.864e + 01 
+ 3.139e + 01 
+ 2.414e + 01 
+ 1.689e + 01 
+ 9.633e + 00 
+ 2.380e + 00 


Fig.3.24 : Contraintes de Von Mises pour d r = 60% (trou horizontal). 



E, Max. Principal 
(Avg: 75o/o) 

+ 9.291e-01 
+ 8.517e-0i 
+ 7.743e-01 
+ 6.968e-01 
+ 6.194e-01 
+ 5.420e-01 
+ 4.645e-01 
+ 3.871e-01 
+ 3.097e-01 
+ 2.322e-01 
+ 1.548e-01 
+ 7.737e-02 
-6.354e-05 


Fig.3.25 : Déformation pour d r = 60% (trou horizental). 
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Type du trou 

Contraint maximale o vm 
(MPa) 

Deformation maximale £ max 

Vertical 

60.14 

0.5 

Horizontal 

89.42 

0.93 

Ecart (%) 

32.74 

46.24 


Tab.3.8: comparaison des valeurs maximales pour les deux structures (trouée 
horizontalement et verticalement). 


Concernant la disposition du trou à la verticale, on peut remarquer que la répartition des 
contraintes dans la structure trouée (figure 3.22) est identique à celle observée dans la 
structure non trouée (figure 3.19), ainsi que les valeurs maximales de la contrainte a vm = 61.2 
MPa pour la structure non trouée contre a vm = 60.14MPa pour la structure trouée. Les figures 
(3.20 et 3.23) montrent que la présence du trou n’influe pas sur la valeur maximale de la 
déformation, elle atteint 50% tandis que les déformations sont beaucoup plus importantes prés 
du trou. 

Concernant la disposition à l’horizontale, il est évident qu’elle influence la direction du trou 
(tableau 3.8). Les contraintes et les déformations sont plus faibles dans l’alignement du trou 
positionné à la verticale ; d’après ces résultats, nous constatons qu’il existe un écart de 
32.74% de contrainte et 46.24% de déformation. Les figures (3.24) et (3.25) montrent que la 
zone voisinant le trou est bien une zone de concentration de contraintes et de déformations. 

3.7. Conclusion. 

Les essais de compression uniaxiale réalisée avec le logiciel des éléments finis ABAQUS ont 
permis de révéler des courbes liées au comportement mécanique du polyuréthane cellulaire ; 
ces courbes présentent toutes la même allure, montée en charge élastique, phase plateau et 
densification, ce qui est conforme à la littérature. La réponse mécanique des solides étudiés et 
en parfait accord avec la théorie de Gibson et Ashby. Notre solution numérique représentée 
par les courbes contraintes déformation et les valeurs de la limite élastique a* t coincide avec 
la solution analytique ; on trouve une différence moins de 0.94 au niveau de la limite 
élastique. Soit un écart inférieur à 4.49% entre la solution analytique et la solution numérique. 
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Les résultats ont révélé que la densité influe sur les valeurs de la limite élastique a* t , la 
déformation de la densification Sd ainsi que, les valeurs maximales de la déformation et la 
contrainte que peut subir le polyuréthane cellulaire. 

Les résultats obtenus montrent que, pour une même intensité de charge, l’énergie absorbée est 
plus importante dans le polyuréthane cellulaire que dans celui massif. Ce qui explique le 
choix des solides cellulaires dans l’emballage et les applications de conforts pour les sièges. 

La présence d’ouverture dans les structures cellulaire et l’orientation des trous influent 
considérablement sur la rigidité du matériau. 
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Conclusion 


Conclusion 

L’étude numérique entreprise dans cette thèse nous a permis de simuler le comportement 
élastoplastique d’un solide cellulaire de polyuréthane à cellules ouvertes en 3D. 

Cette étude décrit avec précision le comportement en compression statique des mousses en 
polyuréthane. 

La transition du matériau de l’état élastique à l’état plastique est gouvernée par l’effondrement 
élastique pour les élastomères, l’effondrement plastique pour les mousses rigides, et par le 
mécanisme de rupture dans le cas des mousses fragiles. 

Un code de calcul basé sur la méthode des éléments finis a été utilisé pour l’analyse de 
structure polymère en élastoplasticité. Sur la base des résultats obtenus qui coïncident bien 
avec les valeurs analytiques (travaux de Gibson et Ashby), on peut conclure que la méthode 
basée sur le modèle Hyperfoam pour l’évaluation des contraintes s’avère satisfaisante. 

Les résultats obtenus mettent en évidence l’importance de la densité, 

La présence d’ouverture dans la structure influe considérablement sur la rigidité du matériau. 

Pour une éventuelle étude complémentaire dans ce contexte, on donne les suggestions 
suivantes : 

-Tenir compte d’autres caractéristiques des mousses comme la morphologie des cellules, leur 
taille, leur arrangement.... 

-Prendre en considération le phénomène d’anisotropie dans le matériau. 

-considérer l’effet du vieillissement et de la température. 

-compléter cette étude par la prise en compte d’autre types de structures, d’étendre l’analyse 
en expérimentale. 
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Résumé 

Le travail présenté dans cette thèse concerne l’étude numérique du comportement mécanique 
des polymères cellulaires, le matériau considéré dans ce travail est le polyuréthane à cellules 
ouvertes. La méthode des éléments finis a été employée pour déterminer les contraintes et les 
déformations dans le matériau. Le calcul consiste plus particulièrement à présenter les 
résultats d’un modèle macroscopique en trois dimensions par éléments finis. On utilise le 
code de calcul par éléments finis Abaqus qui permet de modéliser les mousses polymères par 
le modèle Hyperfoam. 

Mots clés : polymère cellulaire, polyuréthane cellulaire, éléments finis, Abaqus, propriétés mécaniques, le model 
hyperfoam 


Abstract 

The study presented in this thesis relates to the numerical study of the mechanical behavior of 
foam polymers, the material considered in this work is the polyuréthane with open cell. 

The finite element analysis was used to predict stress and deformations in polymer foam 
materials. 

The objective of this Work is to présent results of a macroscopie model in three dimensions 
by finite éléments. 

The hyperfoam model was used for predicting the behavior of polymer foam; this model is 
available in ABAQUS program. 

Keywords: polymer foam, polyuréthane foam, finite element, Abaqus, Mechanical properties, hyperfoam model 
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